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Abstract. An investigation of the stability of the transfer function of the European Southern Observatory's Very Large 
Telescope Interferometer has been carried out through observations of Fomalhaut, which was observed over a range in hour 
angle from 21:50-05:24 on 20 October 2002. No significant variation in the transfer function was found for the zenith angle 
range 5°-70°. The projected baseline varied between 139.7 m and 49.8 m during the observations and, as an integral part of 
the determination of the transfer function, a new accurate limb-darkened angular diameter for Fomalhaut of 2.109 ± 0.013 
mas has been established. This has led to improved values for the emergent flux = (3.43 ± 0. 10)xl0 8 Wm -2 , effective 
temperature = 8819 ± 67 K and radius = (1.213 ± 0.01 l)x 10 9 m (R/R = 1.744 ± 0.016). The luminosity has been found 
to be (6.34 ± 0.20)x 10 27 W (L/Lq = 16.5 ± 0.5). 
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1. Introduction 

A major problem in ground-based optical/infrared long- 
baseline interferometry is establishing the true interference 
fringe visibility from measurements made through the Earth's 
turbulent atmosphere. The observed visibility is always less 
than the true visibility due to residual seeing and instrumen- 
tal effects no matter how well the interferometer is designed, 
compensated and aligned. The technique generally adopted to 
correct for the loss in visibility is to interleave observations 
of the object being studied with observations of a calibrator - 
a star for which the true visibility at the baseline and wave- 
length being employed is known. A calibrator is thus a star 
that is either essentially unresolved or whose angular diame- 
ter is known with sufficient accuracy to predict the true visi- 
bility with the desired accuracy. Division of the observed vis- 
ibility for the programme object by the observed visibility for 
the calibrator yields the true visibility for the programme ob- 
ject, providing the instrumental and atmospheric conditions 
do not change significantly between the measurements. 

Optical/infrared interferometers generally measure the 
square of the fringe visibility as defined by Michelson and, 
for convenience, we will follow Davis et al. ( 1999b I and use 
correlation C to represent the square of the fringe visibility 
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V. We define the transfer function Tp of an interferometer as 
the ratio of the observed correlation C b s over the expected 
correlation C cxp . The transfer function is thus given by 
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The true correlation for a programme object is then given 
by dividing its observed correlation by the transfer function 
for its associated calibrator. 

This procedure works well when programme object and 
calibrator are close in the sky and the observations are made 
close in time. However, problems have been experienced at 
blue wavelengths (^440 nm) with the Sydney University 
Stellar Interferometer (SUSI) (Davis et al. 1999a) when the 
program object and calibrator are separated by more than 
w 5° on the sky, or when the seeing conditions are chang- 
ing on the timescale of the individual measurements. 

As part of the commissioning programme of the Euro- 
pean Southern Observatory's Very Large Telescope Interfer- 
ometer (VLTI) (Glindemann et al. 2000 1, a large number of 
calibrator observations have been accumulated (Percheron et 
al. 2003 ) and are being pursued to assess, among other things, 
the stability and repeatability of the transfer function in the 
spectral K band at a wavelength of 2.2 /xm. The correlation 
measurements are made with the VINCI instrument (Kervella 
et al.|2000 1, which is based on the FLUOR instrument (Coude 
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du Foresto et al. 1998 1, and it employs an optical fibre beam- 
combiner and photometric monitoring of the beams. The spa- 
tial and temporal scales of the turbulence-induced wavefront 
distortions scale as A 6 / 5 and, with the spatial filtering pro- 
vided by VINCI, the VLTI is expected to be significantly less 
affected by atmospheric effects than instruments working at 
shorter wavelengths and with conventional beam-combiners. 

One step proposed for the investigative programme was 
to observe a single star over a large range in hour angle, and 
hence in zenith angle Z, to investigate whether the transfer 
function had a dependence on Z. Time was allocated for these 
programmes when the VLTI was operating with its two test 
siderostats (Derie et al. 2000 1 separated by 140 m. The lim- 
iting K magnitude of the system for an unresolved star was 
~ +3.0 and, for a partially resolved star, the limiting magni- 
tude was brighter by 2.5 logC CX p- With this sensitivity limit a 
shorter baseline would have provided a choice of unresolved 
stars for the experiment. However, given the circumstances, it 
was decided to make observations at the longer baseline since 
they were of intrinsic interest and, in addition, would provide 
experience for planning more detailed experiments with a 
shorter baseline in the future. The limiting magnitude placed 
severe restrictions on the choice of stars, particularly as an 
unresolved star with declination close to the latitude of the 
VLTI location on Paranal (« —25°) would ideally be chosen 
for following over a large range in hour angle. At the sched- 
uled time for the observations no suitable unresolved candi- 
date stars were available. It was decided to follow Fomal- 
haut (a PsA, HR 8728, HD216956, spectral type A3V), even 
though it would be significantly resolved at 140 m based on 
the Narrabri Stellar Intensity Interferometer limb-darkened 
angular diameter of (2.10 ± 0.14) mas (Hanbury Brown et 
al. [T974l . 



2. The Observations and Data Analysis 

The observations were made on the night of 20 October 2002 
in the normal observing mode of the VLTI as a series of Ob- 
serving Blocks (OBs) except that, instead of alternating OBs 
between a programme star and a calibrator, all OBs were for 
Fomalhaut. Each block consisted of three groups of 500 scans 
of the fringe envelope. The star was followed over the hour 
angle range from 21:50 to 05:24 giving a range in zenith an- 
gle from 5° to > 70°. The sky was clear throughout the obser- 
vations and the mean zenith seeing, measured in the visual, 
was 0.76". The seeing only exceeded l"for the first 3 of the 
34 accepted sets of scans. The atmospheric coherence time, 
again determined for the visual, had 99% of the 502 values 
between 1.1 ms and 2.7 ms with a mean value of 1.85 ms. 
The range of values for the seeing and coherence time rep- 
resent a typical night at Paranal and in no sense are they ex- 
ceptional and, as discussed in Sect. |4] there is no correlation 
between the variations in the seeing or coherence time and in 
the observed correlation values. 

The observational data were processed in the VINCI 
pipeline (Ballesteret al. 2002). The processing automatically 
excludes scans for which the signal/noise in the photometry 
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Fig. 1. The correlation measurements as a function of hour 
angle on the night of 20 October 2002. The broken line is the 
expected correlation for a uniform-disk angular diameter of 
2.086 mas (see Sect.|4j. Further details are given in the text. 



channels is less than 5. Some 25% of the sets of scans re- 
tained only a small fraction of the total of 500 scans due to 
drifting of the starlight off the input end of the optical fibers. 
The groups retaining less than 40% of the total scans (15 out 
of a total of 49 sets) were found to deviate randomly from 
the general trend of the observations and to have large uncer- 
tainties associated with them. They were removed from the 
analysis reported here, although a separate analysis showed 
that they had no significant influence on the results or con- 
clusions because of their low statistical weight. 

The correlation values from VINCI's two interferometry 
channels, A and B, differ by a scale factor but track each other 
closely. Initially they were analysed independently and gave 
essentially identical results with no significant differences. 
The two sets of correlation values were therefore put on the 
same scale by normalising the B values by the weighted mean 
ratio of the A/B values. The two sets of data were then com- 
bined by taking a weighted mean of the A and scaled B val- 
ues at each hour angle. The weighted mean values of correla- 
tion as a function of hour angle, zenith angle and projected 
baseline are listed in Table \l\ Figure ^ shows the correla- 
tion values plotted against hour angle for the combined data. 
The large variation of correlation with hour angle was ex- 
pected due to the change in projected baseline from 139.7 m 
to 49.8 m during the course of the observations. 

3. The Effective Wavelength 

In interferometric measurements, the spatial frequency of 
each observation is given by the ratio of the baseline and the 
wavelength. While the baseline can be reliably measured and 
its projection easily predicted as a function of observation 
time, an accurate determination of the wavelength requires 
some effort since the measurements are usually performed 
using filters with a finite bandwidth. The spectral sensitiv- 
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Table 1. The weighted mean values of correlation from the 
two VLTI signal channels as a function of hour angle, zenith 
angle and projected baseline. Note that the formal uncertain- 
ties in the correlation values have been increased by a factor 
of 2.2 as explained in Section^ 
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ity curve of the VLTI VINCI system has been determined 
from the spectral transmission curve of the K band filter used, 
the attenuation curve for the optical fibers, the detector spec- 
tral quantum efficiency curve, and representative transmis- 
sion curves for the atmosphere across the band (Manduca & 
Bell [1972 Lord fT992l . The spectral reflectivities of the sil- 
ver and gold coated mirrors in the optical path are essentially 
constant across the K band. The spectral sensitivity curves 
were computed for the two atmospheric transmission curves 
and gave effective wavelengths, for a constant photon flux 
per unit wavelength across the spectrum, of 2.199 /im and 
2. 195 /im. The mean of these two values, 2. 197 /im, has been 
adopted. Based on the differing model predictions for the at- 
mospheric transmission across the band, the range of precip- 
itable water vapour at Paranal, and the effect of zenith angle, 
the uncertainty in the effective wavelength is estimated to be 
±0.005 /im. More details on this approach can be found in 



a technical report by Davis & Richichi (2003J. For the ob- 
servations of Fomalhaut, taking into account the slope of the 
stellar spectrum across the spectral sensitivity curve, the ef- 
fective wavelength was determined to be 2.181 ± 0.005 /im. 
This value is insensitive to the temperature of the star and 
changing the temperature by ± 2000 K changes the effective 
wavelength by <0.001 /im. 



4. The Transfer Function 

The variation in correlation seen in Fig. ^is primarily due 
to the variation in projected baseline during the course of the 
observations and follows closely the curve expected for a par- 
tially resolved star. If there is a variation in the transfer func- 
tion Tp with time this will also contribute to the variation of 
correlation but to a very much smaller degree. 

Initially Tp was computed as a function of zenith angle 
with the expected correlation C cxp computed using the Inten- 
sity Interferometer angular diameter. The limb-darkened an- 
gular diameter of 2. 10 mas (Hanbury Brown et al.. 119741 was 
converted to the equivalent uniform-disk angular diameter for 
the K band for this purpose using a correction factor (p) taken 
from Davis et al. (2000 1 who computed correction factors for 
Kurucz model atmospheres. For an A3V star, assuming an 
effective temperature of 8750 K and log g of 4.0, the ratio 
of the limb-darkened to uniform-disk angular diameter (p) is 
1.01 1 giving an equivalent uniform-disk angular diameter of 
2.08 mas for the K band. The computed values of Tp were 
almost constant with zenith angle showing only a very small 
linear slope. However, the distribution of the observational 
points about the line of best fit had a greater spread than ex- 
pected for the formal uncertainties in the individual values. 
It was clear that the uncertainties were underestimated and, 
although the cause cannot be identified from the data, short 
term effects due to the atmosphere are a likely significant con- 
tributor. The formal uncertainties were increased by a scaling 
factor and \ 2 was computed for a linear fit to the Tp versus 
zenith angle data. This was repeated for a range of scaling 
factors and it was found that a scaling factor of 2.2 gave a 
X 2 probability of ^70%. This factor has been adopted and 
applied to the formal uncertainties. The scaled uncertainties 
have been used throughout the following analysis and discus- 
sion. It is important to stress that this has no effect on the 
results or conclusions of this paper apart from increasing the 
uncertainties in the derived parameters. 

Although Tp is almost constant with zenith angle based 
on the angular diameter determined with the NSII, a variation 
of Tp with zenith angle cannot be excluded. If a variation is 
present it will affect the angular diameter determined from 
the values of correlation versus projected baseline. Assum- 
ing for the moment that Tp is constant with zenith angle, and 
following the general practice for stars with compact atmo- 
spheres, the values of correlation as a function of projected 
baseline were fitted with the curve for a uniform disk 
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Fig. 2. Normalised observed correlation values versus pro- 
jected baseline. The curve is the fit for the equivalent 
uniform-disk angular diameter of 2.086 mas. 

where x = tvWub /A, b is the projected baseline of the inter- 
ferometer, #ud is the equivalent uniform-disk angular diam- 
eter of the star and A is the effective wavelength of the obser- 
vations. The projected baselines were computed in the VLTI 
data-processing pipeline and the effective wavelength for the 
observations was 2.181 ± 0.005 fim (Sect. [3}- The fit yielded 
a uniform-disk angular diameter of 2.086± 0.012 mas and a 
correlation value at zero baseline of 0.3904 ± 0.0048 with 
a reduced \ 2 value of 0.83. The uncertainty in the angular 
diameter is the formal error of the fit. 

When the uncertainty in the effective wavelength (Sect. [3 
is taken into account, the uniform-disk angular diameter be- 
comes 2.086 ± 0.013 mas. The observed correlation values, 
normalised by the zero-baseline correlation value, and the fit- 
ted equivalent uniform disk curve are shown in Fig.|2] 

It can be seen from Fig.[2]that the fit is excellent with no 
obvious trends away from the uniform disk curve. The points 
at the shortest baselines, which show the greatest scatter and 
largest uncertainties, were taken at zenith angles exceeding 
65°but have negligible effect on the resulting uniform-disk 
angular diameter. There is no correlation between the differ- 
ences of the observed correlation from the fitted curve and 
the variations in the seeing and atmospheric coherence time. 

The transfer function has been assumed constant in the 
analysis. If it had varied significantly during the course of the 
observations it would have affected the correlation values and 
led to an erroneous value for the angular diameter. Assume 
for the moment that Tp varied with zenith angle and that the 
variation is known and used to 'correct' the observed val- 
ues of correlation. An angular diameter fit to the 'corrected' 
values of correlation would be expected to result in an im- 
proved fit with a lower value for the reduced \ 2 - To explore 
this possibility various linear dependencies of the Tp on both 
the zenith angle Z and on sec -3 / 5 Z were tried. The latter 
is the theoretical dependence of the atmospheric spatial and 
temporal coherence parameters ro and to on Z. In each case 
the correlation values were 'corrected' and an angular diam- 
eter fit was made. It was found that the minimum value for 



the reduced \ 2 occurred for zero dependence on Z and on 
sec -3 ^ 5 Z. It is concluded that any dependence of Tp on Z 
is small and that the angular diameter determined on the as- 
sumption of zero dependence is valid. 

Finally, Tp was computed as a function of Z, using the 
new uniform-disk angular diameter to predict the expected 
correlation in equation^ The resulting values for Tp are plot- 
ted against zenith angle in Fig. [3] The figure shows no appar- 
ent dependence of transfer function on zenith angle over the 
range from 5°-70°, as confirmed through a linear regression 
fit. Again there is no correlation between the differences of 
the observed correlation from the fitted curve and the vari- 
ations in the seeing and atmospheric coherence time. The 
mean transfer function was 0.3901 ± 0.0026, in agreement 
with the zero-baseline correlation value obtained from the an- 
gular diameter fit. 

It is noted that in determining the angular diameter, cir- 
cular symmetry has been assumed. The value of V sin i of 
85 kms -1 for Fomalhaut (Abt & Morrell. [T995l l does not clas- 
sify it as a rapid rotator, and any effect of rotational distortion 
on the observations is small compared with the uncertainty in 
the angular diameter. Fig.|4]shows the u,v track of the obser- 
vations. The angle between the projected baseline and a fixed 
arbitrary reference axis through the centre of the star in the 
plane of the sky changed by only 45.8°during the observa- 
tions. Fomalhaut is a main-sequence IR excess star but shows 
no excess radiation at the wavelength of observation (Gillett 
1986 1. Scattered radiation from the star will be a small frac- 
tion of the direct stellar radiation and will not significantly 
affect the results discussed here. An independent study based 
on VLTI measurements of Fomalhaut, using the standard ap- 
proach of calibrated observations, also confirms that there is 
no detectable infrared excess at the VINCI wavelengths (Di 
Folco et al. l2004l 

The mean transfer function corresponds to a fringe vis- 
ibility of 62.5% for a completely unresolved source and it 
only applies to the particular night of observation. The trans- 
fer function may vary due to changes in the atmosphere and 
in the instrument. Long term changes occur for the VLTI due 
to a temperature-dependence of polarization in the fiber beam 
combiner and changes in the optical train to configure differ- 
ent baselines. An investigation of the long-term changes in 
the transfer function of the VLTI is being pursued and val- 
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Fig. 3. The transfer function Tp as a function of zenith angle. 
The line is a linear regression fit to the data. Further details 
are given in the text. 
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Fig. 4. The u, v plot of the observations showing the 40°range 
in orientation of the projected baseline relative to the star. 
Further details are given in the text. 

ues of the fringe visibility between 40% and 87% for an un- 
resolved source have been measured over a period of about 
18 months of operation of the VLTI. The implication of the 
present result is that, in the short term and with a single base- 
line, the VLTI is remarkably stable. 

5. Fundamental Quantities for Fomalhaut 

It was concluded in the previous section that any dependence 
of the Tp on Z is small and that the uniform disk angular di- 
ameter determined on the assumption of zero dependence is 
valid. Based on this assumption, we have proceeded to derive 
the limb-darkened angular diameter and some fundamental 
parameters for Fomalhaut. The factor p, for converting the 
uniform disk angular diameter to the true limb-darkened an- 
gular diameter, has been computed to be 1.011 ± 0.001 for 
the relatively broad spectral passband of the VLTI follow- 
ing the procedure given by Tango & Davis (2002i. Thus the 
true limb-darkened angular diameter is 2.109 ± 0.013 mas. 
This is in excellent agreement with the limb-darkened angu- 
lar diameter determined with the Narrabri Stellar Intensity 
Interferometer of 2.10 ± 0.14 mas and with a factor of >10 
improvement in accuracy from ~ ±6.7% to ~ ±0.6%. 

The combination of the new accurate limb-darkened an- 
gular diameter with the flux received from the star, and 
with the Hipparcos parallax, enables improved values for the 
emergent flux at the stellar surface, and the star's effective 
temperature and radius to be calculated. 

The flux received from Fomalhaut was derived by Code 
et al. J1976I . Following the same procedures, the integrated 
flux has been revised by putting the visual spectrophotome- 
try of Davis & Webb J1974l > on the Vega calibration of Hayes 
fT985t in place of that of Oke & Schild fi970l . and using the 
absolute flux calibration at 550 nm of Megessier ( 1995 ). The 
infrared contribution has been estimated from the flux distri- 



bution for Fomalhaut given by Alekseeva et al. ( 1997 ), scaled 
to match the Davis & Webb flux distribution in the range 700- 



800 nm, the JHKL photometry of Bouchet et al. ( 199 lj con- 
verted to fluxes using Megessier's calibration dl995> . and the 
12 /.im flux given by Gillett ( 1986 1. The ultraviolet contribu- 
tion in the band 182-330 nm has been reduced by 2.5% in 
accord with the change in the visual flux between 330 nm 
and 370 nm to which the uv spectrophotometry was origi- 
nally matched by Code et al. The flux shortward of 182 nm 
has not been changed. The total flux received from the star, 
integrated over the entire spectrum, is estimated to be (8.96 
±0.25)xl0~ 9 Wm- 2 . 

Combining the received flux with the limb-darkened an- 
gular diameter gives the emergent surface flux to be (3.43 
± 0.10)xl0 8 Wm" 2 and the effective temperature 8819 ± 
67 K. The accuracy of the effective temperature has been im- 
proved by a factor of five compared with the determination 
by Code et al. ( 1976 1 and is now limited by the uncertainty 
in the integrated flux rather than by that of the angular diam- 
eter. The improved accuracy will assist in identifying stellar 
atmosphere models with the star. 

The Hipparcos parallax ( fTW7l of 130.08 ± 0.92 mas, 
combined with the limb-darkened angular diameter, gives the 
radius of Fomalhaut as (1.213 ± 0.011)xl0 9 m (R/Rq = 
1.744 ± 0.016). Previous estimates of the radius have been 
based on the Narrabri Stellar Intensity Interferometer angular 
diameter which limited the accuracy to ~ ±6.7% but, with 
the new value for the angular diameter, the accuracy has been 
improved to ±0.9%. The luminosity of Fomalhaut, given by 
47r<i 2 /, where d is the distance and / is the total flux received 
from the star, is (6.34 ± 0.20) x 10 27 W (L/L Q = 16.5 ± 0.5). 

6. Summary 

We have analyzed a set of VLTI interferometric measure- 
ments obtained on the bright star Fomalhaut using the VLTI. 
These measurements were recorded with the main goal of 
studying the accuracy and stability of the VLTI equipped with 
the VINCI instrument. Unlike standard interferometric prac- 
tice, they were not interspersed with independent calibrator 
observations. 

For the night concerned, while it is not possible to com- 
pletely rule out a dependence on zenith angle, we have shown 
that the transfer function was remarkably stable even though 
the observations were taken over large ranges in zenith an- 
gle, projected baseline, and observed correlation values. The 
observational data are an excellent fit to the transform for a 
uniform disk and our analysis shows that, if there is a depen- 
dence on zenith angle, it is at a level that does not significantly 
affect the value of the fitted uniform disk angular diameter. 

A bonus, and an integral result of the analysis of the ob- 
servations, was the accurate determination of the angular di- 
ameter for Fomalhaut. This has allowed fundamental quanti- 
ties for this important star to be determined with significantly 
improved accuracy. 

The angular diameter has been determined without re- 
course to observations of a calibrator. This has been possi- 
ble as a result of observations having been made over a wide 
range in projected baseline length combined with the stability 
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of the transfer function. It is emphasised that the interleav- 
ing of observations of calibrators with those of programme 
stars will continue to be the normal mode of operation, since 
it would not be realistic to abandon it on the basis of one 
night's observations. However, the implications of the results 
presented here are that the VLTI in its current configuration, 
using the VINCI optical fiber beam combiner, gives remark- 
ably stable transfer functions in the short term with a single 
baseline, and that the angular separation of target and cali- 
brators is not as crucial a matter as found at shorter wave- 
lengths with conventional beam-combiners. It augurs well for 
the future of the VLTI and for optical/infrared interferometry 
in general. 
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